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Abstract 

The Rosenblatt process is a self-similar non-Gaussian process which lives in second 
Wiener chaos, and occurs as the limit of correlated random sequences in so-called "non- 
central limit theorems" . It shares the same covariance as fractional Brownian motion. We 
study the asymptotic distribution of the quadratic variations of the Rosenblatt process 
based on long filters, including filters based on high-order finite-difference and wavelet- 
based schemes. We find exact formulas for the limiting distributions, which we then use to 
devise strongly consistent estimators of the self-similarity parameter H . Unlike the case of 
fractional Brownian motion, no matter now high the filter orders are, the estimators are 
never asymptotically normal, converging instead in the mean square to the observed value 
of the Rosenblatt process at time 1. 
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1 Introduction 



Self-similar stochastic processes are of practical interest in various applications, including 
econometrics, internet traffic, and hydrology. These are processes X = {X (t) : t > 0} whose 
dependence on the time parameter t is self-similar, in the sense that there exists a (self- 
similarity) parameter H 6 (0, 1) such that for any constant c > 0, {X (cb) : t > 0} and 
{c H X{t) : t > 0} have the same distribution. These processes are often endowed with other 
distinctive properties. 

The fractional Brownian motion (fBm) is the usual candidate to model phenomena in 
which the selfsimilarity property can be observed from the empirical data. This fBm B H is 
the continuous centered Gaussian process with covariance function 

R H (t, s) := E [B H (t) B H (s)] = l -{t 2H + s 2H - \t - s\ 2H ). (1) 

The parameter H characterizes all the important properties of the process. In addition to 
being self-similar with parameter H, which is evident from the covariance function, fBm has 
correlated increments: in fact, from ([1]) we get, as n —* oo, 

E [(B H (n) - B H (1)) B H (1)] =H(2H- 1) n 2H ~ 2 + o {n 2H ~ 2 ) ; (2) 

when H < 1/2, the increments are negatively correlated and the correlation decays more slowly 
than quadratically; when H > 1/2, the increments are positively correlated and the correlation 
decays so slowly that they are not summable, a situation which is commonly known as the 
long memory property. The covariance structure ([1]) also implies 



E 



B H (t) - B H {s)f =\t-s\ 2H ; (3) 



this property shows that the increments of fBm are stationary and self-similar; its immediate 
consequence for higher moments can be used, via the so-called Kolmogorov continuity criterion, 
to imply that B has paths which are almost-surely (H — e)-H61der-continuous for any e > 0. 

It turns out that fBm is the only continuous Gaussian process which is selfsimilar with 
stationary increments. However, there are many more stochastic processes which, except for 
the Gaussian character, share all the other properties above for H > 1/2 (i.e. ([1]) which implies 
([2]), the long- memory property, ([3j), and in many cases the Holder-continuity). In some models 
the Gaussian assumption may be implausible and in this case one needs to use a different 
self-similar process with stationary increments to model the phenomenon. Natural candidates 
are the Hermite processes: these non-Gaussian stochastic processes appear as limits in the 
so-called Non-Central Limit Theorem (see [5], [8], [25]) and do indeed have all the properties 
listed above. While fBm can be expressed as a Wiener integral with respect to the standard 
Wiener process, i.e. the integral of a deterministic kernel w.r.t. a standard Brownian motion, 
the Hermite process of order q > 2 is a qth iterated integral of a deterministic function with 
q variables with respect to a standard Brownian motion. When q = 2, the Hermite process is 
called the Rosenblatt process. This stochastic process typically appears as a limiting model in 
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various applications such as unit the root testing problem (see [3T] ), semiparametric approach 
to hypothesis test (see [13]), or long-range dependence estimation (see |15j). On the other hand, 
since it is non-Gaussian and self-similar with stationary increments, the Rosenblatt process can 
also be an input in models where self-similarity is observed in empirical data which appears 
to be non-Gaussian. The need of non-Gaussian self-similar processes in practice (for example 
in hydrology) is mentioned in the paper [26] based on the study of stochastic modeling for 
river-flow time series in [16J. Recent interest in the Rosenblatt and other Hermite processes, 
due in part to their non-Gaussian character, and in part for their independent mathematical 
value, is evidenced by the following references: [4], [6], [10], p], [19], [20], [27], [28]. 

The results in these articles, and in the previous references on the non-central limit 
theorem, have one point in common: of all the Hermite processes, the most important one in 
terms of limit theorem, apart from fBm, is the Rosenblatt process. As such, it should be the 
first non-Gaussian self-similar process for which to develop a full statistical estimation theory. 
This is one motivation for writing this article. 

Since the Hurst parameter H, thus called in reference to the hydrologist who discovered 



its original practical importance (see |14j). characterizes all the important properties of a Her- 

mite process, its proper statistical estimation is of the utmost importance. Several statistics 
have been introduced to this end in the case of fBm, such as variograms, maximum likelihood 
estimators, or spectral methods, k-variations and wavelets. Information on these various ap- 
proaches, apart from wavelets, for fBm and other long-memory processes, can be found in the 
book of Beran [3]. More details about the wavelet-based approach can be found in [2], [Tl] and 



In this article, we will concentrate on one of the more popular methods to estimate 
H: the study of power variations; it is particularly well-adapted to the non-Gaussian Hermite 
processes, because explicit calculations can be performed via Wiener chaos analysis. In its 
simplest form, the kth power variation statistic of a process {X t : t G [0, 1]}, calculated using 
iV data points, is defined as following quantity (the absolute value of the increment may be 
used in the definition for non-even powers): 



There exists a direct connection between the behavior of the variations and the conver- 
gence of an estimator for the selfsimilarity order based on these variations (see [7] , [28] ) : if the 
renormalized variation satisfies a central limit theorem then so does the estimator, a desirable 
fact for statistical purposes. 

The recent paper [28] studies the quadratic variation of the Rosenblatt process Z (the 
V/v above with k = 2), exhibiting the following facts: the normalized sequence 
satisfies a non-central limit theorem, it converges in the mean square to the Rosenblatt random 
variable Z (1) (value of the process Z at time 1); from this, we can construct an estimator for 
H whose behavior is still non-normal. The same result is also obtained in the case of the 




[30]. 




k 
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estimators based on the wavelet coefficients (see [2]). In the simpler case of fBm, this situation 
still occurs when H > 3/4 (see for instance |29J). For statistical applications, a situation in 
which asymptotic normality holds might be preferable. To achieve this for fBm, it has been 
known for some time that one may use "longer filters" (that means, replacing the increments 
Xi+i — X_i_ by the second-order increments Xi+i — 2Xj_ + Xi-i , or higher order increments 

N N N N N 

for instance; see [7])- To have asymptotic normality in the case of the Rosenblatt process, it 
was shown in [28] that one may perform a compensation of the non-normal component of the 
quadratic variation. In fact, this is possible only in the case of the Rosenblatt process; it is not 
possible for higher-order Hermite processes, and is not possible for fBm with H > 3/4 [recall 
that the case of fBm with H < 3/4 does not require any compensation]. The compensation 
technique for the Rosenblatt process yields asymptotic variances which are difficult to calculate 
and may be very high. 

The question then arises to find out whether using longer filters for the Rosenblatt 
process might yield asymptotically normal estimators, and/or might result in low asymptotic 
variances. In this article, using recent results on limit theorems for multiple stochastic integrals 
based on the Malliavin calculus (see [22] , [23] ) , we will see that the answer to the first question 
is negative, while the answer to the second question is affirmative. We will use quadratic 
variations (k = 2) for simplicity. A summary of our results is as follows. Here denotes the 
underlying probability space, and L 1 (fi) and L 2 (Q) are the usual spaces of integrable and 
square-integrable random variables. 

• Vn = Th + T4 where is in the ith Wiener chaos (Proposition [2]) . 

• ^~^4 converges in distribution to a standard normal (Theorem [2]), where c\^jj is given 
in Proposition HI 

• N , Vn and N , T2 both converge in L?(Q) to the Rosenblatt random variable Z(l) 

(Theorem [3]); the asymptotic variance C2,h is given explicitly in formula (I16p in Propo- 
sition 

• There exists a strongly consistent estimator Hn for H based on Vn (Theorem [5]), and 

2 c 2 ]^ 2 (log N)N l ~^ N (^Hn — H^j converges in L 1 (£1) to a Rosenblatt random variable 

(Theorem [7|). Here c%n is again given in (|16p . Note that while the rate of convergence 
of the estimator, of order N~ 1+H log -1 N, depends on H , the convergence result above 
can be used without knowledge of H since one may plug in Hn instead of H in the 
convergence rate. 

• The asymptotic variance c^h i n the above convergence decreases as the length of the 
filter increases; this decrease is much faster for wavelets-based filters than for finite- 
difference-based filters: for values of H < 0.95, C2,h reaches values below 5% for wavelet 
filters of length less than 6, but for finite-difference filters of length no less than 16. 

/C2.H 



When H G (1/2,2/3), then y N — Z(l) converges in distribution to a stan- 
dard normal, where c^h is given explicitly in formula (116j) and C3 u i n formula (I19p . 
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Similarly, for the estimator we have that 



-2log(H N -H)-^Z(l) 



converges 

in distribution to the same standard normal. However, no mater how much we increase 
the order and/or the length of the filter, we cannot improve the threshold of 2/3 for H. 

What prevents the normalization of Vm from converging to a Gaussian, no matter how 
long the filter is, is the distinction between the two terms T2 and T4. In the case of fractional 
Brownian motion, Vn contains only one "T2"-type term (second chaos), but this term has a 
behavior similar to our term T4, and does converge to a normal when the filter is long enough; 
this fact has been noted before (see [7]). In our case, the normalized T2 always converges (in 
I? (O)) to a Rosenblatt random variable; the piece that sometimes has normal asymptotics 
is T4, but since T2 always dominates it, Vjv's behavior is always that of T2. This sort of 
phenomenon was already noted in [6] with the order-one filter for all non-Gaussian Hermite 
processes, but now we know it occurs also for the simplest Hermite process that is not fBm, 
for filters of all orders. 

The organization of our paper is as follows. Section 2 summarizes the stochastic analytic 
tools we will use, and gives the definitions of the Rosenblatt process and the filter variations. 
Therein we also establish a specific representation of the 2-power variation as the sum of two 
terms, one in the second Wiener chaos, which we call T2, and another, T4, in the fourth Wiener 
chaos. Section 3 establishes the correct normalizing factors for the variations, by computing 
second moments, showing in particular that T2 is the dominant term. Section 4 proves that 
the renormalized T4 is asymptotically normal. Section 5 proves that T2 converges in I? (ft) 
to the value Z (X) of the Rosenblatt process at time 1. In Section 6 it is shown that the 
variation obtained by subtracting this observed limit of T2 leads to a correction term which is 
asymptotically normal. Section 7 establishes the strong consistency of the estimator H for H 
based on the variations, and proves that the renormalized estimator converges to a Rosenblatt 
random variable in L 1 (O). Its asymptotic variance is given explicitly for any filter, thanks to 
the calculations in Section 3. In Section 8, we compare the numerical values of the asymptotic 
variances for various choices of filters, including finite-difference filters and wavelet-based filters, 
concluding that the latter are more efficient. 



2 Preliminaries 

2.1 Basic tools on multiple Wiener-Ito integrals 

Let {Wt : t e[0, 1]} be a classical Wiener process on a standard Wiener space (0, J 7 , P). If a 
symmetric function / G L 2 ([0, l] n ) is given, the multiple Wiener-Ito integral I n (f) of / with 
respect to W is constructed and studied in detail in [211 Chapter 1]. Here we collect the results 
we will need. For the most part, the results in this subsection will be used in the technical 
portions of our proofs, which are in the Appendix. One can construct the multiple integral 
starting from simple functions of the form / := ^ c^...^!^ x ...xA in where the coefficient 
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c h,..,i„ is zero if two indices are equal and the sets Ai j are disjoint intervals, by setting 

Uf):= E c ll ,..,„^(^ 1 )...W(^J 

it,..., in 

where we put W {\ a ,b]) = W([a, 6]) = W> ~~ Wo 5 then the integral is extended to all sym- 
metric functions in L 2 ([0, l] n ) by a density argument. It is also convenient to note that this 
construction coincides with the iterated Ito stochastic integral 

In(f)=n\ / .../ f(t u ...,t n )dW tl ...dW tn . 
Jo Jo Jo 

The application I n is extended to non-symmetric functions / via 

In{f)=In{!) (5) 

where / denotes the symmetrization of / defined by f(x\ , . . . , x x ) = ^ YlaES n f( x v(i) x a(n) ) 
The map (n!) _1 ^ 2 / n can then be seen to be an isometry from L 2 ([0, l] n ) to L 2 (S7). The nth 
Wiener chaos is the set of all integrals {/„(/) : / G L 2 ([0, l] n )}; the Wiener chaoses form 
orthogonal sets in L 2 (O). Summarizing, we have 

E(I n (/)Im(ff)) ="!(/, 5)i2([0,i]«) ifm = n, (6) 
E(I n {f)I m {g)) = iim + n. 

The product for two multiple integrals can be expanded explicitly (see |21|): if / G 
L 2 ([0, l] n ) and g G L 2 ([0, l] m ) are symmetric, then it holds that 

mAn 

Uf)l m (g) = £ £!C4^J m+n _ 2 ,(/ ® / 5 ) (7) 

where the contraction f (E>e9 belongs to L 2 ([0, i] m + n ^ 2 ^) for £ = 0, 1, ... ,m An and is given by 

e, h, . . . , t 

m—t J 

= / f(si,...,s n -.e,ui,...,u e )g(ti,...,t rn ^£,u 1 ,...,u e )dui...du£. 

J[0,i\e 

Note that the contraction (/ <g)£ g) is not necessary symmetric. We will denote by (f®ig) its 
symmetrization. 

Our analysis will be based on the following result, due to Nualart and Peccati (see 
Theorem 1 in [22J). 

Proposition 1 Let n be a fixed integer. Let I n (/7v) be a sequence of symmetric square inte- 
grable random variables in the nth Wiener chaos such that limN-*ooE [l n (fN) 2 ] = 1- Then 
the following are equivalent: 
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(i) As N — ► oo ; the sequence {I n (fN) '■ N > 1} converges in distribution to a standard 

Gaussian random variable. 

(ii) For every r = 1, . . . , n — 1 

2.2 Rosenblatt process and filters: definitions, notation, and chaos repre- 
sentation 

The Rosenblatt process is the (non-Gaussian) Hermite process of order 2 with Hurst index 
H G (|, 1). It is self-similar with stationary increments, lives in the second Wiener chaos and 
can be represented as a double Wiener- ltd integral of the form 

zW{t) := Z(t) = f f L t ( yi ,y 2 )dW yi dW y2 . (8) 
Jo Jo 

Here {Wt,t G [0, 1]} is a standard Brownian motion and L t (yi,y2) is the kernel of the Rosen- 
blatt process 

where 

tr = S±l and ah) = -±- ( , H X 1 ' 2 

2 y ' H + l \2{2H-\) J 

and K H is the standard kernel of fBm, defined for s < t and H G (^, 1) by 

K H (t,s):=c H s^- H J (u- s) H ~lu H -^du (10) 
i 

where ch = ( p(2-2H~H-±) ) ^ an< ^ ') * s ^ e ^ e * a f unc ti° n - F° r * > s, we have the following 
expression for the derivative of K H with respect to its first variable: 

^(t, s) := d 1 K H (t, s) = c H (£) (t - . (11) 

The term Rosenblatt random variable denotes any random variable which has the same distri- 
bution as Z(l). Note that this distribution depends on H. 

Definition 1 A filter a of length £ G N and order p G N \ is an (Z + 1)- dimensional vector 
a = {ao, Oil, . . . ,af\ such that 

i 

^ a q q r = 0, /or 0<r<p-l, rGZ 

I 

9=0 
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with the convention 0° 



If we associate such a filter a with the Rosenblatt process we get the filtered process 
V a according to the following scheme: 



V c 



: =zZ a i Z 



q=0 



N 



for i = £,..., N - 1. 



Some examples are the following: 
1. For a = {1,-1} 



ya 



N 



N 



This is a filter of length 1 and order 1. 
2. For a = {1,-2,1} 



z — 

N 



2Z 



N 



+ Z 



i - 2 
N 



This is a filter of length 2 and order 2. 



3. More generally, longer filters produced by finite-differencing are such that the coefficients 
of the filter a are the binomial coefficients with alternating signs. Therefore, borrowing 
the notation V from time series analysis, VZ (i/N) = Z (i/N) — Z ((i — 1) /N), we define 
V-? = VV J_1 and we may write the jth-order finite-difference-filtered process as follows 

V» (%) := (V'Z) (±) . 



N 



From now on we assume the filter order is strictly greater than 1 (p > 2). 

For such a filter a the quadratic variation statistic is defined as 



V N : = 



1 

N-e 



N-l 

E 



" \v a (i?)\ 



1 
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Using the definition of the filter, we can compute the covariance of the filtered process V a (^) : 

'i+f 



N 



a q a r E 



q,r=0 
N -2H 



Z 



N 



N 



Z 



i + j — r 
N 



q,r=0 
N -2H ' 



a q a r (\i - q\ 2H + r\ 2H -\j + q- r\ 2H ) 

N -2H ' 



^2 Oi q a r \j + q 



r\ 2H + 



a q a r {\i-q\ 2H + \i+j-r\ 2H ). 



q,r=0 q,r=0 

Since the term r=0 a q a r (\i — q\ 2H + \i + j — r\ 2H ^ vanishes we get that 

N -2H ' 



^2 a q a r \j + q 



\2H 



(12) 



q,r=0 

Therefore, we can rewrite the variation statistic as follows 



V, 



1 



N 



N -£ 

2N 2H 
N-l 

2N 



JV-l 

E 



vr-(O) 



1 



c(H)(N-l) 



^ a r a q \q - r 

q,r=0 
■ill N-l 



2H 



i=e 



N-l 

E 



*&(0) 



N 



where 



1 2H 



c(H) = — ^2 ®rOt q \q — r\ 

q,r=0 

The next lemma is informative, and will be useful in the sequel. 



(13) 



Lemma 1 c(H) is positive for all H G (0, 1]. Also, c(0) = 0. 

Proof. For H < 1, we may rewrite c(H) by using the representation of the function 
\q — r\ 2H via fBm B , as its canonical metric given in ([3]), and its covariance function Rjj 
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given in (pTJ) . Indeed we have 

I 

c(x)=-J2 a raq B\(B H (q)-B H (r)Y 



q,r=0 



/I a r oi q {Rh (q, q) + Rh (r, r) - 2R H (q, r)) 



q,r=0 



-2 y^aJ E a r -R# (r, r) +2 E a r a q Rn (q, r) 



,9=0 



vr=0 



<j,r=0 



+ 2 E a r a q R H (q,r) = E 

q,r=0 



E (?) 

>9=0 



> 



where in the second-to- last line we used the filter property which implies X^o=o a i = 0> anc ^ 
the last inequality follows from the fact that Y2 q =o a qB H (q) is Gaussian and non-constant. 
When H = 1, the same argument as above holds because the Gaussian process X such that 
X (0) = and E (X (t) - X (s)) 2 = \t - s\ 2 is evidently equal in law to X (i) = tN where N 

is a fixed standard normal r.v. The assertion that c(0) = comes from the filter property. ■ 
Observe that we can write the filtered process as an integral belonging to the second 
Wiener chaos 



q=0 



N 



h ( ^2 a q L igL ) : = h (Ci) , 

,q=0 



where 



(14) 



Ci '■— / ttqLi-q. 
L — ' N 

q=0 

Using the product formula ([7|) for multiple stochastic integrals now results in the Wiener chaos 
expansion of Vn ■ 

Proposition 2 With Ci as in jl4\ ), the variation statistic Vn is given by 



2N 



2H N-l 



N 



c(H)(N-l) 



E 



2N 



2B 



c{H){N-l) 
:=T 4 + T 2 , 



7V-1 



N-l 



E h {Ci ® d) + 4 E h {Ci ®i Ci 



i=£ 



where T4 is a term belonging to the 4th Wiener chaos and T 2 a term living in the 2nd Wiener 
chaos. 
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In order to prove that a variation statistic has a normal limit we may use the charac- 
terization of JV(0, 1) by Nualart and Ortiz-Latorre (Proposition [T]) . Thus, we need to start by 
calculating E [| V}v 1 2 ] so that we can then scale appropriately, in an attempt to apply the said 
proposition. 

3 Scale constants for T2 and T4 

In order to determine the convergence of Vn, using the orthogonality of the integrals belonging 
in different chaoses, we will study each term separately. This section begins by calculating the 
second moments of T2 and T4. 

In this section we use an alternative expression for the filtered process. More specifically, 
denoting b q := Ylr=o a r> we rewrite Cj as follows, for any % = £,..., N — 1: 



Ci t i := Ci = y^^QgL 



q=Q 



N 



N N 



a o ( Li — Li-i ) + («o + «i) ( Li-i — Li-2 )+... + («o + . . . + a^_i) ( — L 



N N 



q=0 



N 



L i — q 



(15) 



Recall that the filter properties imply Ylq=o a g = an d act = — J2q=o a q- 
3.1 Term T 2 

By Proposition O we can express E(T 2 2 ) as: 



2 _ 64 N 4H 

E(T2) " c(H)HN-ey 



E 



'JV-1 



£ h {Ci ®i Ci) 



2! 64 N 4H ^ 
c(H) 2 {N — I) 2 ^ \ ' 01 " j 01 j 'L 2 ([o,i] 2 ) 



Proposition 3 We have 



lim E 

jV— »oo 



\N 1 - H T 2 f 



C2,H- 



where 



C2,H 



64 



2H-1 



c(H) 2 \H(H + 1) 2 



bqb r 



q,r=0 



\l + q-r\ 2H + \l-q + r\ 2H - 2\q - r\ 2H 



(16) 



This proposition is proved in the Appendix. 
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3.2 Term T 4 



In this paragraph we estimate the second moment of T4, the fourth chaos term appearing in the 
decomposition of the variation Vjy. Here the function ^dS^ 1 (C, ® Cj) is no longer symmetric 
and we need to symmetrize this kernel to calculate T^s second moment. In other words, by 
Proposition [21 we have that 



E (Tf ) 



AH 



c(H) 2 (N - if 
c(H) 2 (N-£) 2 



E 



'N-1 



N-1 



([o,i] 4 ) 



where Ci®Ci := Cj <S> C{. Thus, we can use the following combinatorial formula: 
If / and g are two symmetric functions in L 2 ([0, l] 2 ), then 

4!(/®/,5®5}z,2([o,i]4) 

= (2!) 2 (/ ® f,g ® <7>L2([o,i]4) + (2!) 2 (/ 0i 5,9 <8>i /^([o.lp)- 



It implies 



E (T|) 



4/\T 



JV-l 



r4! X)(C , i®C , i J C7 ; -®C ; -) ia ([o I i]*) 

AT-l 

c(g) 2 (iV-^) 2 4 E^^^i^ilW) 
N-1 



c{H) 2 {N -if 
4N iH 



+ 



c(H) 2 (N-£)' 



([o,i] 2 ) 



:= T 4) ( 1) + 14,(2)- 



The proof of the next proposition, in the Appendix, shows that the two terms Z^i) and Z4,(2) 
have the same order of magnitude, with only the normalizing constant being different. 

Proposition 4 Recall the constant c (H) defined in lilty) . Let 



oo £ „ 

n,H ■= Yl b qi b q2 b ri b r2 j 



dudvdudv' 



k=i qi,q2,ri,ri=0 

\u — v + k — qi+ r\ 
\u - v! + k - qi + q 2 



2H'-2 I / I . , . 

\U — V + k — q2 + T2 



\2H'-2 



2H'-2 



v — v + k — r\ + r 2 



2H'-2 
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and 



p%{k) :-- 



\2H 



c(H) 



Then we have the following asymptotic variance for y NT4 : 

2 



lim E 

AT— >CX> 



NTa 



c hH :=4! [l + Y,\p a H( k )\ 2 )+Ti,H- 



(17) 



fc=0 



This proposition is proved in the Appendix. Observe that in the Wiener chaos decom- 
position of V/v the leading term is the term in the second Wiener chaos (i.e. T2) since it is of 
order N 11 ^ 1 , while T4 is of the smaller order iV -1 / 2 . We note that, in contrast to the case of 
filters of lenght 1 and power 1, the barrier H = 3/4 does not appear anymore in the estimation 
of the magnitude of T4 Thus, the asymptotic behavior of V/v is determined by the behavior of 
T2. In other words, the previous three propositions imply the following. 



Theorem 1 For all H £ (1/2, 1) we have that 



lim E 

N^oo 



N 



l-H 



v N 



C2,H, 



where C2,h is defined in |7 



^From the practical point of view, one only needs to compute the constant C2,h to 
find the first order asymptotics of V/v- This constant is easily computed exactly from its 
formula (|16p . unlike the constant c± t H in Proposition (jj]) which can only be approximated via 
its unwieldy series-integral representation given therein. 



4 Normality of the term T4 

We study in this section the limit of the renormalized term T4 which lives in the fourth 
Wiener chaos and appears in the expression of the variation Vn- Of course, due to Theorem 
[T] above, this term does not affect the first order behavior of Vjy but it is interesting from the 
mathematical point of view because its limit is similar to those of the variation based on the 
fractional Brownian motion (|29j). In addition, in Section 6, we will show that the asymptotics 
of T4, and indeed the value of ci t H, are not purely academic. They are needed in order to 
calculate the asymptotic variance of the adjusted variations, those which have a normal limit 
when H £ (1/2,2/3). 

Define the quantity 

N„ \fN 2N 2H 



Gn := —T A = y J2 T * ^ ® C ^ 

ci,h ^/c^H c{H){N - I) ^ 



Tt 9 N 2H W-l \ 

■ - T X;(Ci®Ci) :=h(jBN)- (18) 



c(H) (N 



i=£ 
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^From the calculations above we proved that lirnjv-»ac E(G^) = 1. Using the Nualart-Peccati 
criterion in Proposition Q3 we can now prove that Gat is asymptotically standard normal. 

Theorem 2 For all H £ (1/2, 1) Gn defined in J7^|) converges in distribution to the standard 
normal. 



Setup of proof of Theorem [2J To prove this theorem, by Proposition H] and Propo- 
sition [H it is sufficient to show that for all r = 1,2,3, 

JV^So W 9N ® t9N Wl 2 ([0,1]( 8 - 2 ^) = °- 

For t = 1, 2, 3, this quantity can be written as 



lim 



< lim 



lim 



4N 4H+1 


c hH c(H) 2 (N- 








ci,hc(H) 2 (N — 


I) 2 


4N m+1 




c liH c(H) 2 {N- 


I) 2 



N-l 



^(Cj<8> Ci)<§ T (Cj :® Cj] 



N-l 



(d ® Ci) ® T (Cj ® Cj 



1,3=1 
2 N-l 



L 2 ([0,1]( 8 ~ 2t )) 
2 



L 2 ([0,1]( 8 - 2 ^)) 



((Ci ® Ci) ® T (Cj ® Cj), (C m ® C m ) ® T (C n ® C n )) . 



i,j,m,n=£ 



The Appendix can now be consulted for proof that for each r = 1,2,3 this quantity 
converges to 0, establishing the theorem. □ 



5 Anormality of the T2 term and Asymptotic Distribution of 
the 2-Variation 

For the asymptotic distribution of the variation statistic we have the following proposition. 

Theorem 3 For all H G (1/2,1), both To and the normalized 2-variation —. Vn con- 
verge in L 2 (Q) to the Rosenblatt random variable Z(l). 

Setup of proof of Theorem [3]. The strategy for proving this theorem is simple. First 
of all Proposition [4] implies immediately that ./V 1_ ^T4 converges to zero in L 2 (Q). Thus if we 
can show the theorem's statement about T2, the statement about Vn will following immediately 
from Proposition [2j 

Next, to show ^ C2 H ffz converges to the random variable Z (1) in L 2 (0), recall that T2 is 
a second-chaos random variable of the form ^(In), where /iv(yi, 1/2) is a symmetric function in 
L 2 ([0, l] 2 ), and that this double Wiener-Ito integral is with respect to the Brownian motion W 
used to define Z (1), i.e. that Z (1) = I2 (L\) where L\ is the kernel of the Rosenblatt process 
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at time 1, as defined in Q. Therefore, by the isometry property of Wiener- ltd integrals (see 
([6])), it is necessary and sufficient to show that ^==f^ converges in L 2 ([0, l] 2 ) to L\. This is 
proved in the Appendix. □ 



6 Normality of the adjusted variations 

In the previous section we proved that the distribution of the variation statistic V/v is never 
normal, irrespective of the order of the filter. However, in the decomposition of Vn, there 
is a normal part, T4, which implies that if we subtract T 2 from Vjv the remaining part will 
converge to a normal law. But T 2 is not observed in practice. Following the idea of the adjusted 
variations in [28], instead of T 2 we subtract Z(l) which is observed. Z(l) is the value of the 
Rosenblatt process at time 1. Thus, we study the convergence of the adjusted variation: 

yN-^Z(l) = V N -T 2 + T 2 -^Z(l) 
:= T 4 + U 2 . 

In Section 4 we showed that -^-T^ converges to a normal law. For the quantity U 2 we 
prove the following proposition 



Proposition 5 For H G (^f), \/NXJ 2 converges in distribution to normal with mean zero 
and variance given by 

c 3 ,h ■= c 2 , H - k - l)k 2H F (£) , (19) 

where c 2y n is defined as in \16\) and F is defined as follows 

e r 

F(x) = d(H) 2 a(H) 2 ^ / dudvdu'dv'\(u-u +q 2 - qi )x + l\ 2H '~ 2 

128a(H) 2 d(H) 2 , , ]2H ,_ 2 . , , t , 2H '-2 



c 2 , H c( H ) 2 



\u — v — q\ + r\\ \u — v — q 2 + r 2 \ 



\(v - v - r\ + r 2 )x + 1 7— . -\u — v — qi + n 

y/C2JlC(H) 

\(v - u - q 2 + ri)x + l\ 2H 2 + \(u - u + qi - q 2 )x + l\ 2H 2 

Proof. The proof follows the proof of |28|, Proposition 5] and is omitted here. 
Therefore, for the adjusted variation we can prove the following 
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Theorem 4 Let Z t : t £ (0, 1) be a Rosenblatt process with H £ (1/2, 2/3). TTien i/ie adjusted 
variation 

-^tea)-^!)). 
converges to a standard normal law. Here c^h, c^Hi and c^^h are given in |77| ), [TS\) , and 

Proof. The proof follows the steps of the proof of [251 Theorem 6] and is omitted. ■ 



7 Estimators for the self-similarity index 

We construct estimators for the self-similarity index of a Rosenblatt process Z based on the 
discrete observations at times 0, jj, jr, . . . , 1. Their strong consistency and asymptotic distri- 
bution will be consequences of the theorems above. 



7.1 Setup of the estimation problem 

Consider the quadratic variation statistic for a filter a of order p based on the observations of 
our Rosenblatt process Z: 

i=e y g =o v 7 J 

We have already established that E [Sn] = — ^2 — Yfqr=o a q a r\q — f\ 2H (see expression (fT2j) 
). By considering that E [Sn] can be estimated by the empirical value Sn, we can construct 
an estimator Hn for H by solving the following equation: 

N -2H N * 

Sn = — 2^ a i a r\ a ~ r \ 2HN ■ 

q,r=0 

In this case, unlike the case of a filter of length 1 which was studied in |28j, we cannot com- 
pute an analytical expression for the estimator. Nonetheless, the estimator Hn can be easily 
computed numerically by solving the following non-linear equation for fixed N, with unknown 
xe [1/2,1]: 

J\f-2x 1 

— Yl «9«rk - r\ 2x - S N {2, a) = 0. (21) 

q,r=0 

This equation is not entirely trivial, in the sense that one must determine whether it has a 
solution in [1/2,1], and whether this solution is unique. As it turns out, the answer to both 
questions is affirmative for large N, as seen in the next Proposition, proved further below. 

Proposition 6 Almost surely, for large N, equation [2l\) has exactly one solution in [1/2,1]. 
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Definition 2 We define the estimator Hn of H to be the unique solution of $21\). 

Note that Equation (|2ip can be rewritten as Sn = c(x)N~ 2x /2 where the function c 
was defined in (|13p . The proposition is established via the following lemma. 

Lemma 2 For any H £ (1/2, 1), almost surely, Huin^oo N 2H Sn = c(H) /2. 

Proof. Firstly, we show that Vn converges to zero almost surely as N —* oo. We 
already know that this is true in L 2 (O). Consider the following 

P(\V N \ > N~P) < N qf3 E(\V N \ q ) < c qA [E(V$)] q/2 <cN ql3 N^ H -^ q . 

If we choose (3 < 1 — H and q large enough so that (1 — H — (3)q > 1. This implies that 

oo oo 

£ P (VjvI > iV^) < c £ N^+ H - 1 ^ < +oo 

AT=0 iV=0 



Therefore, the Borel-Cantelli lemma implies |Vjv| — > a.s., with speed of convergence equal to 
N~ p , for all /3 < 1 - H. Since V N = - 1 we have 

2N 2H 

1 + V N = -—i : — S N = 2N 2H S N /c (H) . (22) 

E q , r =o a q a r \q-r\ 2H 

The almost-sure convergence of Vn to yields the statement of the lemma. ■ 

Proof of Proposition [6l For x € [^,1] and for any fixed N, define the function 

Fn (x) = ^N~ 2x -S N = Yl a ^\l ~ r \ 2X ~ s *- 

q,r=0 

Equation (f2Tj) is Fn (x) = 0. Observe that Fn{x) is strictly decreasing. Indeed, we have that 



F' N (x) = log (N- 2x ) a?Orl9 " r ? X ~ N' 2x ^ a q a r log \q - r\ 

q,r=0 q,r=0 

Then, F' N (x) < is equivalent to 

r Eg r=0 a Q a r lo § I? ~ A \l ~ r ? X 1 



r\ 2x . 



N > exp 



E ? ,r=o <V*r|<Z-r 



2x 



since we know, using Lemma[U that c (x) = r=0 a g a r |g — r\ 2x , which is evidently continuous 
on [5, 1], is strictly negative on that interval. Thus, if we choose N to be large enough, i.e. 



N > max exp 



11 II l9r 

L g ,r=o a q a r log |q - r\ \q - r\ zx 



^ih 1 ] [ T, q ,r=o a q a r\q- r 



q,'i 
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the function Fn (x) is invertible on [|,1], and equation (|21|) has no more than one solution 
there. 

To guarantee existence of a solution, we use Lemma[2j This lemma implies the existence 
of a sequence sn such that 2N 2H Sn = c{H) + en and \miN^oo£N = almost surely. Since in 
addition c is continuous, then almost surely, we can choose N large enoug h, so that 2N 2H S N 
is in the image of [|, 1] by the function c. Thus the equation c(x) = 2N 2H Sn has at least 
one solution in [i, 1]. Since this equation is equivalent to (|21|) . the proof of the proposition is 
complete. ■ 

7.2 Properties of the estimator 

Now, it remains to prove that any such Hn is consistent and to determine its asymptotic 
distribution. 

Theorem 5 For H G (1/2, 1) assume that the observed process used in the previous definition 
is a Rosenblatt process with Hurst parameter H . Then strong consistency holds for Hn, 

lim Hn = H, a.s. 

AT— >oo 

In fact, we have more precisely that lirriTv^oo ( H — Hn ) log N = a.s. 



Proof. ^,From line (|22p in the proof of Lemma [21 and using the fact that H n solves 
equation (HQ, i.e. c (h n ) N~ 6n = 2S N , we can write 



i , T/ _ 2N<2H c c ( H n) M 2(H~H N ) 

Now note that c{Hn)/c(H) is the ratio of two values of the continuous function c 
at two points in [1/2,1]. However, Lemma Q] proves that on this interval, the function c 
is strictly positive; since it is continuous, it is bounded above and away from 0. Let a = 
niiaggM / 2) i] c(x) > and A = max^g^^i] c(x) < oo. These constants a and A are of course 
non random. Therefore c(Hn)/c(H) is always in the interval [a/A, A/a]. Thus, almost surely, 



log [c(H N )/c(H] 



, A 
< log—. 

a 



log(l + VN) = 2lH-H N )logN + log(^^-) . (23) 



We may now write 



c{H) ^ 

Since in addition lirriAr^oo log (1 + Vn) = a.s., we get that almost surely, 



H-H 



N 



o 



1 

logiV 



18 



This implies the first statement of the proposition. 

The second statement, which is more precise, is now obtained as follows. Since Hn — * H 
almost surely, and c is continuous, log [c{Hn)/c{H)^ converges to 0. The second statement 
now follows immediately. ■ 

The asymptotic distribution of the estimator Hn is stated in the next result. Its proof 
uses Theorem [3] and Theorem [IJ plus the expression (|23p . While novel and interesting, this 
proof is more technical than the proofs of the proposition and theorem above, and is therefore 
relegated to the Appendix. 

Theorem 6 For any H G (i,l), the convergence 

lim 2c~]l 2 N 1 - H (h n - h) logiV = Z(l) 

holds in L 2 (Q), where Z(l) is a Rosenblatt random variable. 

As can be seen from Theorem [3] and Theorem [UJ the renormalization of the statistic 
Vjv, as well as the renormalization of the difference Hn — H, depend on H: it is of order 
of N^ H . The quantities iV 1_ - ff Vjv and N l ~ H Hn cannot be computed numerically from the 
empirical data, thereby compromising the use of the asymptotic distributions for statistical 
purposes such as model validation. Therefore one would like to have other quantities with 
known asymptotic distribution which can be calculated using only the data. The next theorem 
addresses this issue by showing that one can replace H by H n in the term N l ~ H , and still 
obtain a convergence as in Theorem El this time in L 1 (O). Its proof is in the Appendix. 

Theorem 7 For any H G (i, 1), with the Rosenblatt random variable Z (1), 



lim E 



2 c 2 JfiV 1 -^ log n(h n -h)-Z (1) 



8 Numerical Computation of the Asymptotic Variance 

In practice certain issues may occur when we compute the asymptotic variance. The most 
crucial question is what order of filter we should choose. Indeed, from (|16|) with Hn instead of 
H, it follows that the constants of the variance not only depend on the filter length/order (£, 
p), but also on the number of observations (N). We measure the "accuracy" of the estimator 
Hn by its standard error which is the following quantity: 

2ATi-^iviogAr' 

There are several types of filters that we can use. In this paper, we choose to work with 
finite-difference and wavelet- type filters. 
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• The finite-difference niters are produced by finite-differencing the process. In this case 
the filter length is the same as the order of the filter. The coefficients of the order-^ finite 
difference filter are given by 



• The wavelet filters we are using are the Daubechies filters with /c-vanishing moments. 
(By vanishing moments we mean that all moments of the wavelet filter are zero up to a 
power). The Daubechies wavelets form a family of orthonormal wavelets with compact 
support and the maximum number of vanishing moments. In this scenario, the number 
of vanishing moments determines the order of the filter and the filter length is twice the 
order. For more details, the reader can refer to |17j . 

We computed the standard error for N = 10, 000 observations, filters of order varying 
from 2 to 20 and Hurst parameters varying from 0.55 to 0.95. This means that the correspond- 
ing lengths of the finite-difference filters were 2 to 20 and for the wavelets 4 to 40. The code 
we use to simulate the Rosenblatt process is based on a Donsker-type limit theorem and was 
provided to us by J.M. Bardet pp. The results are illustrated in the figures 1, 2, and 3, on 
the next page; these are graphs of the asymptotic standard error ^/C2 : h / (2N 1 ~ Hn log N) for 
various fixed values of H as the order of the filters increase. 

We observe that the standard error decreases with the order of the filter. Furthermore, 
we observe that the wavelet filters are more effective than the finite-difference ones, since they 
have a higher impact on the decrease of the standard error for the same order, as the filter 
increases. Specifically, the graph in Fig. 1, with the finite difference filters, shows that for fixed 
H, there is no advantage to using a filter beyond a certain order p, since the standard error 
tends to a constant as p — > oo. This does not occur for the wavelet filters, where the standard 
error continues to decrease as p —* oo in all cases as seen in the graph in Fig. 2. On the other 
hand, the finite-difference filters have lower errors than the wavelet filters for low filter lengths; 
only after a certain order p* do the latter become more effective; this comparison is seen in 
the graph in Fig. 3, where p* is roughly 9. 

In addition, since the order of convergence depends on the true value of the Hurst 
parameter H, we investigated the behavior of the error with respect to H. It seems that the 
higher H is, the more we lose in terms of accuracy; this is visible in all three graphs. 

In general, the choice of a longer filter might lead to a smaller error, but at the same 
time it increases the computational time needed in order to compute H and its standard error. 
In a future work, we will study extensively this trade-off and other consequences of using longer 
filters. 




fc+i 
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Standard Error vs. Order of Filter (Finite Difference) 
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Fig. 1: Finite Difference Filters. 
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Fig. 2.: Wavelet Filters. 
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Fig. 3.: Comparison between the two types of filter. 
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9 Appendix: proofs. 

9.1 Proof of Proposition [3] 

We start by computing the contraction term Cj (g>i Cj: 



{Q (8>i Ci) (2/1,2/2) = / Ci(x,yi)Ci(x,y 2 )dx 
Jo 

e .1 

= S~] b q b r / (L i-h-i) (x,yi) - Li-q(x,yi)) IL i-(r-i) (x,y 2 ) - Li-r(x,y 2 ) ) cfc 

g,r=0 

I i-g+l . i-r+1 

9 \ / JV " JV 

= d(H) >^ b g 6 r l [0 ^-g+i ] fai)l [0 »-H-i ] (j!/9) / fife 

« v,_n 



d ( H f ^2 b i brl [o, i ~ q+1 ] ( m ^[o, '-^ 1 ] ^ - 



g,r=0 

dK H ' ( dK H ' dK H ' N 

^ dudv—^-{u,y 1 )—^—(v,y 2 )dudv\J dx — (u,x) — (f,x) 
a{H)d(H) 2 b A l [ 0j *=3±± ] (Vl) 1 [o, i=^i ] (^2) 

rfu du « - v — — (w,yij— — {v,y 2 )dudv, 
ou ov 



where I, 



i— g i— q+1 

g ~~ V "FT' jv 
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Now, the inner product computes as 

(Ci <8>i Ci, Cj <8>i Cj)^2jo,i] 2 

= a{H) 2 d{H) 4 ^2 b qi b ri b q2 b r2 / dy 1 dy 2 
dudvdu dv \u — v\ \u — v \ 



I i qi J U ri J Ijq 2 JI jr 2 

dK H ' dK H ' dK H ' , dK H ' , 

du ( u,yi ' dv ^ V,m ' du ' dv' ^ > y2 > dudvdu dv 

i 

= a(Hfd{H) A Y, KKKK 

qi,ri,q2,T2=0 



/ / / / dudvdu dv \u — v\ 2H ' 
Jii Q1 Jh ri Ji jq2 Ji jr2 

( f uAu ' dK H ' , ^8K H ' , , w \ / r vAv ' dK H> , ^dK H> , , N , ' 
(H) 4 d(H) 4 ^ b qi b ri b q2 b r2 / / / / dudvdu dv 

91,7-1, 92,7*2=0 Jtiqi , ' J *ri ^ 



= a 



w | \2H -2\ ' 1 \2H -2\ ' \2H -2\ ' \2H -2 

X \U — V\ \U — V \u — u\ \V — V \ 



We make the following change of variables 



" 1 " _ '-jr 1 v 
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and the second moment of T 2 becomes 
E (T|) 

N— 1 

b 01 6 ri b q2 b r , 2 / / / dudvdu dv 



128 a(H) 4 d(H) 4 N 4H " 1 



c(H) 2 (N - I) 2 

V ' V 7 i,J=€ <21,ri ,92 ,7-2=0 " Jl <7l " Jl n , ' J Jq 2 "-^'rj 



w | \2H -2\ ' '\2H-2\ ' \2H -2i '|2if-2 

x m-u u — in u — u k> — in 



128 a{H) A d{H) A N 4H 

^j 1 (N-e) 2 ~N4N*H'-s ^ ^ " qi " ri " q2 " r2 J l0 i] 

V 7 V ' i,j=€ ?i,ri,?2,r 2 =0 



AT-1 

dudvdu dv 



/ , / 



X |u- w - qi +n\ 2H 2 \u - v -q 2 + r 2 \ 2H 2 

x |u - u + « - 3 - qi + g 2 | \v — v + 1 — j — ri + r 2 \ 



4j/m4 l n-1 e 

— —2 b qi b ri b q , 2 b r , 2 I dudvdu dv 

' i,j=eqi,ri,q2,r 2 =0 1 "' " 



128 a{Hfd{H) 



x\u-v-qi+n\ 2H 2 \u - v' - q 2 + r 2 \ 2H 2 

x (|u - u + i - j - q 1 + q 2 \ 2H ~ 2 \v - v + i - j - n + r 2 \ 2H . 

Let est. = 128 a ^H)^ H ^ ■ We- study first the diagonal terms of the above double sum 
E { T 2-diag) 



v ' oi .ri .O9.r9=0 l u ' J 



dudvdu dv 

91^1,92,7-2=0 

2Jj'-2i ' ' I |2if'-2i . ' I |2Jf'-2i„ . ' „ , „ |2if'-2 



X |u — v — + ri| \u — v — q 2 + r 2 \ \u — u — q\ + q 2 \ \v — v — r\ + r 2 \ 
We conclude that 

E {Tl diag ) = O (N- 1 ) 
Let's consider now the non-diagonal terms 

t 

E i T 2-off) = 2cst - b qi b ri b q2 b r 2 
qi,ri,q2,T2=0 

\2H' -2 



x / dudvdu dv x \u — v — q\ + ri\ \u — v — q 2 + r 2 \ 
J[0A] 4 



N-l 



(N-e) 2 



\u-u +i-j -q 1 +q 2 \ 2H ~ 2 \v-v' +i-j -n+r 2 \ 2H 2 . J (24) 
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Observe that the term ([23]) can be calculated as follows: 
1 JV-l 

(N _ 2^ \u -u +i - j - qx+n\ \v-v +i- j-ri+r 2 \ 
- N-lN-i 

= TK? — 7a y) X) I" - «' + fc - 91 + 92| 2H '" 2 k - v + k - n + r 2 \ 2H '~ 2 
I N ~ 1 

—jyz J2(N-k-l)\u-u + k-qx + q 2 \ 2H '~ 2 \v -v +k-n + r 2 \ 2H '- 2 



(N-t) 

N-l , 

Z 1 

k=l v 



k=l 

N 



(N - if 



k + l 



N 



u — u k q\ — q 2 



N N N 



2H -2 



v — v k T\ — r 2 



N N N 



2H -2 



We may now use a Riemann sum approximation and the fact that 4H — 4 = 2H — 2 > — 1. 
Since i is fixed and q\ and q 2 are less than £, we get that the term in ([24"]) is asymptotically 
equivalent to 



AT-l 

£ 

k=i 

We conclude that 







2_ff -2 




2H -2 „ 


H) 












iV 




iV 


Jo 



{l-x)x 2H - 2 dx + o(l) 



2H(2H - 1) 



+ o(l). 



2H-2 



x \^ b qi b ri b q2 b r2 / dudvdu dv \u — v — q\ + r\\ 2H ~ 2 \u —v —q 2 + r 2 \ 2H ~ 1 
„_ -Ao.il 4 



gi,n,92,r2=0 
Using the fact that 



[0,1] : 



|n — v — <7 + r| 2 ^ 2 dudv 



1 



2H'(2H' - 1) 
the proposition follows. 



l + <7-r| 2H +|l-g + r| 2iJ -2|q-r 



2H 



|2if 
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9.2 Proof of Proposition [4] 



9.2.1 The term E (T% {1) 
We have 



E (V 2 



jv-i 



' 1 1 ) ~ c(H) 2 {N -I) 2 4 ' ^ ^ ® ° h ° j ® ^'MtM 4 ) 

4AT4H ^-1 



c(H) 2 (N-£) 2 



i,j=£ 



The scalar product computes as 



(Cj,Cj) L 2([ 0) l]2) — ( ^CfcgLi-^, y^arLj-r 



AT * » AT 

\g=0 r=0 



L2([0,l] 2 ) 



ny^a g Li-2 (yi,y 2 ) y^a r L 1 ~r(y 1 ,y 2 ) dy 1 dy 2 
\U N \t^ N J 



d(H) 2 a g a r 

g,r=0 



/7 1 



'-9 

JV 



dtf*' dK H ' 



N 



dK H ' 8K H ' 



Jo 



J/lVj/2 

e 

= d(H) 2 OL q a r 

q,r=0 

e 

= a{H) 2 d{H) 2 a Q a r 

q,r=0 

where a(H) = H{H 2 +1) = H' (2H' - 1) and 

1 



N IN 



dy\dy 2 

uAv q k h' q k h> 

du (u, 2/i ) dv (v, y\)dyi 



1 3— r 
JV /Ar 



JO 



u — t>| 2 ^ 2 dudv 



. '— g J— r 
AT / at 



u — vl 2 ^ 2 dudv 



'0 Jo 



H(2H - 1) 



TV 



2ii 



+ 



3 — r 



N 



2H 



j — i + q — r 



N 



2j 



28 



Using the fact that a ^2H^i) = \ anc ^ the sca l ar product becomes 



I{2H-1) 



ajHf d{H? < 



1 1 

\ E 



q,r=0 



q,r=0 



I — q 



N 



2H 



+ 



N 
j -r 



211 



+ 



J -r 



N 



2H 



j — i + q — r 



N 



2H 



N 



211 



j — i + q — r 



N 



o=0 



N 



2H 



EM+ E 



i 

E< 

q,r=0 



i — j + q — r 



N 



\r=0 
2H " 



\r=0 



J - r 



N 



2H 



2/P 



o=0 



- a 1 a r 



q,r=0 



i — j + q — r 



N 



2H 



The last equality is true since E^„ =0 a q = by the filter definition. Therefore, we have 
JV-l 



Cj)i 2 ([0,l] : 

i,j=e 



E 

N -4H 



q,r=0 
-AH 



i — j + q — r 



N 



211 



N-1N-2 

I EE 

i=e k=o 



E< 

<j,r=0 



k + q — r 



N 



211 



|2ff 



q,r=0 



N-1N-2 

+ iEE 



i=£ fc=l 



E 



q,r=0 



+ g — r 



AT-2 



(H)2 JY-"(JV-<-l) + l ^- (JV _ fc _ 2) 



fc=0 



g,r=0 



k + q — r 



c(H) 



.(N-l- l)N~ 4H N~ 4H+X 



+ 



N~ 



-AH 



N-2 

E 

k=0 



N-2 

E 

k=0 
2 



A r 



i2.fi" 



A r 



2iY 



2_H 



a g a r + 

<j,r=0 



,2H 



+ 



g,r=0 
N -AH 



N-2 

E* 

fc=0 



q,r=0 



\2H 



At this point we need the next lemma to estimate the behavior of the above quantity. This 
lemma is the key point which implies the fact that the longer variation statistics has, in the 
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case when the observed process is the fractional Brownian motion, a Gaussian limit without 
any restriction on H (see |12j). 

• Lemma 3 For all H 6 (0, 1) ; we have that 



(i) EfcLl Eg,r=0 + 9 " 

(") J2T=l k Yfn,r=0 a q a r\ k + ' 



\2H 



\2H 



< + OO 
2 

< +OO. 



Proof. Proof of (i). Let f(x) = r =o a q a r (1 + (l — r)x) 2H , so the summand can be 
written as 



jTa g a r \k + q -r\ 2H = k* H f(l). 

q,r=0 ^ ' 



Using a Taylor expansion at xq = for the function f{x) we get that 



(l + (q- r)xf H m 1 + 2H(q - r)x + . . . + 



2H(2H - 1) . . . (2H - n + 1) 



(g-r) n x n . 



For small x we observe that the function f(x) is asymptotically equivalent to 

2H(2H - 1) . . . (2H - (p - l))x 2p , 

where p is the order of the filter. Therefore, the general term of the series is equivalent 
to 

{2Hf{2H -I) 2 ... (2H - (p - l)) 2 k 4H ~ 4p 

Therefore for all H < p — \ the series converges to a constant depending only on H. Due 
to our choice for the order of the filter p > 2, we obtain the desired result. 

Proof of (ii). Similarly as before, we can write the general term of the series as 



k 



a q a r \k + q — r\ 

q,r=Q 



2H 



k 



k 2H f 



« {2H) 2 (2H - l) 2 . . . (2H - (p - l)) 2 k iH - 4p ~ 1 
Therefore for all H < p the series converges to a constant depending only on H. 
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Combining all the above we have 
4 N AH 



E 



4iV 4H 



c(tf) 2 (iV-f) 2 
N 



4! 



JV-2 



i /v- 4H+1 - 

\c{h?(n-i-i)n-^+^— y: 



k=0 



a^ov | A; + g — r\ 



2H 



q,r=0 



4! 



E 

fc=0 



^ a g a r \ k + q - 

q,r=0 



\2H 



+ 



N 



-AH 



N-2 

E* 



k=0 



a q a r \ k + q — r 



2H 



N 



N l 1 
2- 



(iV-^) 2 V(^-^) 2 (iV-f) 2 



<3,r=0 
JV-2 



E 

fc=0 



a g a r |fc + g — r 

q,r=0 



2H 



N-2 



+ 



+ 



(N-e) 

4! 
^HO 2 

1 



E fe 

fe=0 

C (F) 2 

/V- 

E* 



g,r=0 



2.H 



N 



l + l 



{N-l) 2 (N-£) 2 J (N-lf 



N-2 ^ 

fc=0 



a g a r |/c + q — r 

q,r=Q 



2H 



N-2 



(n - ey 



k=0 



a,a r |fc + (? — r 



2H 



q,r=0 



4! 



c{Hf 



N-2 



c(H) 2 (N- 1 -(£+ 1)N- 2 ) + (N- 1 - 2N~ 2 ) ^ 



fc=0 



y^ a g a r |/c + q — r 

q,r=0 



2H 



N-2 



+ n- 2 Y, k 



k=0 



y~] a g a r |/c + g 

q,r=0 



\2H 



Since the leading term is of order N we have that 



E (r 4 2 {1) ) ~ 4! c(H)- 2 N- 



N-2 



c{h? + y: 



k=0 



y~^ a g a r |/c + g — rj 

g,r=0 



2/f 



If we define the correlation function of the filtered process as 



ir H (k) _ Y,ir=o a <i a r\k + q-r 



2H 



ir H (0) c(H) 

we can express the asymptotic variance limjv^oo N E (t 2 ^-^ in terms of a series involv- 
ing p%{k). 
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9.2.2 The term E \T% {2) 

In order to handle this term we use the alternate expression (|15p of Cj. Therefore, following 
similar calculations as in the T2 case we get that 



E T, 



4,(2); - (A T_£)2 



..r,.r 1= • / [°> 1 



dudvdu ' dv 



9i,92,n,ri=0 



jV-1 

*E 

i,j=i 



2H'-2 I / / 



•u — « + ? — j — q\ + ri \u -v +1- j - q 2 + r 2 



\2H'-2 



\2H'-2 



\u - u + i - j - qi + q 2 \ \v — V + l — J — 7*1 + r 2 \ 



l2i/'-2 



-4,i? 



(JV - £) 2 

N-lN-i-i 

E E 



V b qi b q2 b ri b T2 I 



dudvdu 1 dv 1 



(Jl,(j2,n,ri=0 



\u — v + k — q\ + r\ 



1 2//' 2 

w — it + — gi + ^2 p — v + fc — ^1 + T2 



[0,1H 
2H'-2 I / / 



it — v + k - q 2 + r 2 

2H'-2 



\2H'-2 



7V-£ 



V b qi b q2 b ri b r2 I 

,.r,.r,=0 • / [°> 1 



?i,<j2,n,ri=o 



£(J\T-fc-l) 



£i=0 



dudvdu 1 dv 1 



|n — i> + Jc — gi + n n— t> + & — (72 + 7*2 



|2i/'-2 



|2ff'-2 



u — n+/c — (?i + (?2 w — v + k — r\ + r 2 



2H'-2 
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We study the convergence of the above series as N — > oo 

N-l 



k=0 

\u — u' + k — qi + q 2 

N-l 



\u — v + k — qi + ri\ \u — v + k — q 2 + r 2 



\2H'-2 



2H'-2 



v — v ' + k — r\ + r 2 



2H'-2 



= (^-i)E 



k=0 



2^"' 2 1/ / 

\u — v + k — qi + r±\ \u — v + k — q 2 + r 2 



\2H'-2 



\u — u + k — qi + q 2 

N-l 



2H'-2 



fc=0 



f — V ' + — T\ + T2 



\u — v + k — qi + r\\ \u — v + k — q 2 + r 2 



2H'-2 



\2H'-2 



\2H'-2 



v — v + k — ri + r 2 



2H'-2 



\u — u' + k — qi + q 2 \ 

■■= (i) + m 

Therefore the general term of the series is asymptotically equivalent to 
(2H' -2)...(2H' -2p-l) x4 



(2p)\ 



\2p 



(u - v - qi + ri) p (u - v - q 2 + r 2 ) 



(u - v! - Ql + q 2 f p (v - v' - n + r 2 ) 2 f fc^" 4 "* 



which converges for all H £ 1). We treat the second series (II) in the same way and we get 
that it is asymptotically equivalent to est. k AH ~ 4 ^ 8p . Combining all the above we have 



E 



( T 4,(2)) 



-4, if 



(N-£) 



i 

E KKKK i 



dudvdu'dv' 



gi,92,n,ri=0 
7V-1 



(*-*)E 



[0,1] 4 

2ff'-2 I / / 



|u — u + — q\ + ri | kt — w + fc — q 2 + r 2 



2H'-2 



\u — u + k - qi + q 2 
N-l 



2H'-2 



k=e 



v — v' + k — r\ + r 2 



u — w + /c — u— w+A; — 52 + ^2 



2H'-2 



\2H'-2 



\2H'-2 



\u — u + k — qi + q 2 \ \v — v + k — r\+ r 2 



\2H'-2 
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The leading term in E [T^^-A is of order N 1 and the constant computes as 

oo I „ 

ti,h = z2 b qi b q2 b ri b r2 I 



fc=£ 9i,g2,ri, ri=0 

|u — u + A; — gi + ri 
In - u' + k - q\ + <?2 1 



dudvdu dv 



2H'-2 I / / I t 

ki - u + k - q 2 + r 2 



\v — v' + k — r\ + r2\ 



2H'-2 



Therefore, combining the two terms we get the statement of the proposition. 
9.3 End of proof of Theorem [2] 

Recall that we only need to show that for r = 1, 2, 3 the terms | \gx ® T ^iv| lx, 2 ([o i] 8 ~ 2r ) conver g e 
to as N tends to infinity. 

• Term for r = 1. 



Ji 



c hH c(H) 2 (N - t) 2 

4^4^+1 



2 N-l 



Yj d C i ® Ci) ®1 (Q ® C j)> (C'm ® C m ) ®i (C„ 8) C n )} 



i,j,m,n= 
2 N-l 



c 1)H c(H) 2 (N - ty 

t,j,m,n—<. 

x(Ci 01 Cj,C rn 01 C' n ) i 2Q 0) i]2). 



(Cj, C m ) L 2([ 0) l]2) (Cj, C„,) L 2([ 0i l]5 
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Thus, we have 



Ji 



N-l 



N 8H+2 1 
< CSt.— — — 2^ 



' (n - iy iv 4 



i — m + qi — r\ 



E 



bq± b ri bq 2 b T2 bq z b rz bq 4 b T4 



i,j,m,n=£ qi,ri,q2,r2,qs,r3,q4,,r4,=0 



2H 



N 

dudvdu'dv' 



j -n + q 2 -r 2 



N 



2H 



[0,1]« 



u - v + i - j - q 3 + r 3 



u — u + i — m — q% + 



< est. 



N 



N 

2 N-l 



N 

2H'-2 



2H'-2 



U — V + m — n — (?4 + T& 



V — V + j — n + T3 + T4 



AT 



AT 

2if'-2 n 



2H'-2 



(n - ty 



E E 



&n bq 2 br 2 bq% b r3 bq 4 b r - 4 



i,j,m,n=£ gi,?"l,g2,?"2,g3,»*3,<?4,»*4=0 
\2H I • „ , „ „ \2H 



x |i — m + q\ — r\\ \j — n + q 2 — r 2 \ 



I dudvdv! dv'\u — v + i — j — q% + r^\ 2H 2 \u — v' + m — n — q^ + r^\ 2H 1 

•Ao.ii 4 



L^[o,i] 



x I It — It + 2 — m — 53 + 94 1 |f — f + j — n + T3 + T4| 



As in the computations for T 4j ( 2 ) we can show that the above series converges and thus 
Ji = C(A^- 2 ), which implies that for all H € (§, 1) 

lim Ji = 0. 

N^oo 

Term for r = 2 



j 2 



- r 



VC1,HC(^) 2 (JV-^) 2 
4iV 4i/+l 



2 N-l 

® C) ®2 (Cj ® C,-), (Cm ® C m ) ® 2 (C„ ® C n 

i,j,m,n=£ 
2 N-l 

Cl H c(H) 2 (N -if) ^ ^' ^'^([o.i] 2 ) ( Cm ' C '™)i 2 ([o,i] 2 ) 

' i,j,m,n=£ 
'■{pi, C m ) L 2([ 0)1 ]2) (Cj, C„) L 2([ 0i l]2). 
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N 8H+2 N ~ l 

^ 2 - CSt ' (M _ ¥A X/ Cj}L 2 [0,l] 2 (Cii C m ) L 2yQ^2(C m , C n ) L 2[ 0>1 ]2 (Cj, C n ) L 2[ 0)1 ] 

i,j,m,n=£ 



est. 



(iV - £) 4 



JV-1 



E E 



i,j,m,n=i gi(?2<?3<74=0 



i-j + QI-Q2 



i-m + qi- q 3 



est. 



N 
N 2 



2H 



m- n + q 3 - g 4 



iV 



2fl 



A/ 



j -n + q 2 - g 4 



2ii 



A r 



2H 



JV-1 



X] agiag2 a <73 a <74 I* ~ 3 + 11 - 12 



2H 



(TV - 1) 4 ^ ^ ''/■•■''n 

i ,j,m,n=£ qi(/2<J3<?4=0 



X |i - m + gi - g 3 | 2H \m-n + q 3 - q^\ 2H \j-n + q 2 - qi\ 2H ■ 

The series converges for all H £ (1/2, 1), so the whole term is of order 0(N~ 2 ) which means 
that goes to zero as N — > oo. 



Term for r = 3. 



AN 



AH+X 



c hH c(H) 2 (N-£) 2 



c ltH c(H) 2 (N - I) 2 
x(d ®i Cj,C m 8>1 C n ) 



2 JV-1 

((Ci Ci) 03 (Cj Cj), (C m 8 C m ) 8>3 (C n 8) C„)> 

i ,j,m,n=£ 
2 N-l 

(Cj, Cj) L 2([ 0i1 ]2) (C m ,C n ) i 2([ ,l]2) 

i ,j,m,n=£ 



With similar computations as in the case of T4 we conclude that J3 = 0(A^ 



9.4 Proof of Theorem H 

According to our previous computations we can write 
/jv (2/1,2/2) = 



8AT 



2H AT-1 



c(H){N -I) 



8d(i/) 2 a(^) ^ , „ , x 



c(fT) (N - 1) 



i=£ q,r=0 



[ [ dudv\u-v\ 2H '~ 2 d 1 K H '(u,y 1 )diK H \v,y 2 ) 



36 



Let us show first that we can reduce this function to the interval y\ G [0, % -^\ and y2 G [0, 
We will show that if y\ G Ii q ,y 2 £ [0, '^jf] (and similarly for the situations y\ G [0, ^jf],V2 £ ij r 
and yi G h q -,y2 G ^v) the corresponding terms goes to zero as iV — > oo. We have, due to the 
fact that the intervals l{ q are disjoint, 

N 1-H N 2H N -! * 

II (N _ e) E E Mri/ l9 (yi)i [0 ,^](y2) 

v ; i=e q,r=0 

dudv\u-v\ 2H '- 2 d 1 K H '(u,y 1 )d 1 K H '(v,y 2 )\\ 2 L2m]2} 

N 2+2H N 



= (N- IV E E b qi b ri b q2 b r2 f Iff 



dv'du'dvdu 

%r 2 



x (ju — v\ • \u' — v'\ • \u — u'\ ■ \v — 

N 2+2H | 



2H'-2 



(N -£) 2 N^ A^4(2g'-2) E E MnMr 2 dudvdudv 



x |u - f - + ri| 2H 2 |u' - t/ - 52 + r 2 | 2// 2 
|u - u - gi + q2\ 2H '~ 2 \v - v - n + r 2 | 2// '~ 2 x iV 1 - 2 ^ 

which tends to zero because 2ii~ > 1. 

This proves the following asymptotic equivalence in L 2 ([0, l] 2 ): 

„. 8d(H) 2 a(H) N 2H 

f N (vi,y2) * — — W ^ T) &Al [0 ,V](yi)l[o,^](y2) 

v ; v ' i=t q,r=0 

X 



/ f dudv\u-v\ 2H '~ 2 d 1 K H '(u,y 1 )d 1 K H '(v,y2) 



lq " *i r 



We will show that the above term, normalize by N± H , converges pointwise for y±, y2 G [0, 1] 
to the kernel of the Rosenblatt random variable. 

On the interval Ii q x Ii r we may attemp to replace the evaluation of d\K H at u and v 
by setting u = (i — q)/N and v = (i — r)/N. More precisely, we can write 

d 1 K H '(u,y l )d 1 K H '(v,y2) = Uk h ' (u, yi ) - d 1 K H 'C-j^-, Vl )j d 1 K H '(v,y 2 ) 

+ d 1 K H \^, yi )(d 1 K H \v, y2 )-d 1 K H ' -d 1 K H \ i -^,y 2 ) 

and all the above summand above can be treated in the same manner. For the first one, 
using the definition of the derivative of K H with respect to the first variable, we get for any 
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d 1 K H '(u,y 1 )-d 1 K H 'C- 1 ^-,y l ) 

= c H y{ \{u- yi ) H 2u H 2 - \—- -yij (-^-) 2 

<cuy{ {-JT-K) [ u 2 -(^) 2 

1 / • \ ff-l 

<^ 2 ( ^ 2/1 ) («~(-^)) 2 



2 



and for any y 2 G [0, ^] 

diK h ' {v,y 2 ) = c H y 2 " {v -y 2 ) H ~^v 



3 



< c ^y 2 -jiT - - 2/i ' f — - — ^ 2 ■ 



As a consequence of the above estimates, 



2H at-1 e 



E E Mrl^]^!) 1 ^]^) 



JV-£ 

i=£ g ,r=0 




^Zq " liq 



dvdu\u - v\ 2H '- 2 (dyK H \u, yi ) - dtK^C-^-^yi)] d 1 K H '(v,y 2 ) 



L _ n i+hN-i t 



JV 

i=£ q,r=0 



jv y v ^ / w 
i 7V_1 * 

z=£ g,r=0 



\" 5 A-r \^ 2 /*-r + l. H _i 

yi ( ^ 2 - 



jv y i jv " y v jv 



The quantity ^ 1 [0) i^l)^]^) " yi)^ § " 2/ 2 f ~ § C^ 1 )^ 

parable, for large JV, to the integral fy lVy2 (, u ~ Ui) H ^{ u — V2) H ~ 1 u H ^2 and the term JV 
in front gives the convergence to zero of ([26]) for any fixed yi,y 2 - 
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This means we have proved the following pointwise asymptotically equivalent for /jv(yi, V2)'- 
N^ H 8d(H) 2 a(H) ^ ' 

V^ fNim,m) * y&5<H) (N-i) g W]<*> W](^ 



Recall that 

dvdu\u — v\* AA * = — 

Thus we get 

N l - H 



jf ^ <focfa| M - t>| 2g '~ 2 = 2 ^ (2 ^ 1} {|1 - g + r\ 2H ' + |1 + g - r\ 2H ' - 2\q - rf H '| . 



In (yum) 

£ b q b r {\l- q + r \ 2H ' + \l + q-r\ 2H ' -2\q-r\ 2H '} 



*d(II) 2 «(H) r,, ^ xm 



C2,H c(H) 

' v ' q,r=0 

N-l 



{N-£) 1 v N ,yiJ 1 y N 

Further, we can ignore the terms q/N and r/N in comparison with i/N in the last line above, 
and thus invoke a Riemann sum approximation, which proves that, for every yi,j/2 S (0, l) 2 

n i-h 

= 8d(g)2 ^ } V b q b r {\l-q + r\ 2H ' + \l + q-r\ 2H ' -2\q- r\ 2H '\ 

C2,H C(H) ^ I J 

' v ' q,r=0 

N-l 

= L 1 (y 1 ,y 2 ). 



To finish the proof it suffices to check that iV 1 H /n is a Cauchy sequence in L 2 ([0, l] 2 ). Up 
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to a constant depending on H we have that for all M, N, 

\\ Nl ~ H fN - Ml ~ H fAl\\L 2 ([0,l} 2 ) 



N 2 ~ 2H \\f 

N 2H 



^Nl 2 ([o,i] 2 ) 

+2 t 



+ M 2 2jFf ||/M|||2 ( [ 0il ]2) - 2N 1 H M 1 H (f N , /a/)l 2 ([o,i] 2 ) 



est. 



{N-if 



r; ,„ „ _ ' _n -/I* ,/jJV JjP 



dudvdu'dv' 



est. 



i,j=i <n>n, <?2 ,r 2 =o 
i 

u — ■ 

M 2H+2 M^l 



x\u- v\ 2H '- 2 \u' - v'\ 2H '- 2 \u - u'\ 2H '- 2 \v - V 



91 Tj J?2 Jr 2 

l\2H'-2 



(M - If 



E E 

i,J=^ 91>'*l>?2,» - 2=0 



6 ri 6 92 6 r2 



M J JM JjM 

ii ir i ii2 i r 2 



dudvdu'dv' 



w I \2H'-2\ I /\2H'-2\ i\2H'-2\ /\2H'-2 

X \U — V\ \U — VI \U — U IV — V \ 



est. 



M 1+H N 1+H 

{M -£)(N -£) 



N-lM-l 



KKKK / / / / 



dudvdu'dv' 



i=l j=£ q\,r\,q2,T2=0 



M 

1 " 'Jq 2 ' ~3r 2 



„ I \2H'-2\ I i\2H'-2\ i\2H'-2\ i\2H'-2 

X \u — VI \u — V \ \u — U I IV — V I 



The first two terms have already been studied and will converge to the same constant as 
M, N — > oo. Concerning the inner product, by making the usual change of variable we have 



(MN) 



H+l 



(NM) 



2H'-2 N-lM-l 



(M-£)(N-£) N 2 M 2 



EE E 

i=t j=£ <ji,ri,<j2,»"2=0 

x \u-v - qi+ n\ 2H '~ 2 \u -v' - q 3 + r 3 \ 2H '~ 2 



dudvdu'dv' 



[0,l]4 



N M N N N N 



2H'-2 



v v i j r\ V2 
N~M + N~N~N + N 



2H'-2 



For large i, j we can ignore the terms jSj, j?, ^7, etc., compared to 4r and 4?. Therefore, the 
above quantity is a Riemann sum that converges to the same constant as the squared terms, 
as M, N — > oo. This finishes the proof of the theorem. 

9.5 Proof of Theorem © 

We wish to show that, as N — > oo, 



E := E 



Z(l)-2c$V"*(j 



H n -H)1orN 



0. 



A minor technical difficulty occurs when Vn is not small. We deal with this as follows. We 
decompose the above expectation E according to whether or not \V^\ < 1/2: we have E = 
E\ + Ei where 



E-\ — E 



1/2 AT 1-H 



l\v N \>i/2[Z{l)-2c 2 ^N 



H N -H) IorN 
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Dealing with this term first, Schwarz's and Minkowski's inequalities yields 



E x < 2P 1 / 2 [\V N \> 1/2] ( E 1 



1/2 



Z(l) 



+ Aq l H N 2 - 2H log 2 N E 1 / 2 



H N -H) 



Since Hn is bounded, the sum of the two rooted expectation terms above is bounded above 
by a constant multiple of N 2 ~ 2H . Therefore to deal with E\, one only needs to show that 
P [|Vjv| > 1/2] <C N~ A+AH . It is well known that any random variable X which can be written 
as a finite sum of Wiener chaos terms up to order q satisfies, for any integer n, E [X 2n ] < 
K n ^ q (E [AT 2 ])™ where K n>q depends only on n and q. This can be proved iteratively by using 
formula (|7|), for instance. Therefore, since Vat is a sum of terms in the second and 4th chaos 
(q = 4), by Chebyshev's inequality, and using Theorem [U with N large enough, 



2n 



P[|1V|>1/2]<4"E \\V N \* 
< 8 n K nA cl H N 



< 4 n c n , 4 ( E 

2Hn-2n 



\Vn\' 



It is thus sufficient to choose n = 3 to guarantee that E\ — ► 0. 

We now only need to study Ei- We invoke the mean value theorem to express (^Hn — H 
log AT more explicitly. For any x,y € [1/2, 1], there exists £ E (x,y) such that 

log ^-A- = (x-y) (logc)' (£) . 

c(y) 

Here the function (log c) is bounded on [1/2, 1], because d is bounded and c is bounded below. 
Therefore, denoting by £at G [1/2, 1] the value corresponding to x = H and y = Hn, and using 
line ()23p in the proof of Theorem [5j we can write 

log (1 + V N ) = (h n - H) (2 log N + (log c)' (Cn)) 

log (1 + V N ) 



and thus 



H N -H) (21ogiV) = log(l + -^v) 



21og^ + (logc)'(£ iV ) 

Since |(logc)' (£jv)| is bounded (by a non-random value), by choosing iV large enough, an 
upper bound for the last fraction above, in absolute value, is 2V/v/logA r . Therefore (using 
Minkowski's inequality), 



% = B 1 / 2 
< E i/2 

+ E l/2 



l\v N \<i/2 (Z(l)- 2c 2 ]i 2 N 1 ~ H [H N - H) logiV 



l|^|<i/ 2 ( Z (1) - C2 ' /2 AT 1 - H (log (1 + W)) 



L |Viv| 



< 1/2 (2c 2 j/ 2 iY 1 - H y JV /log 



(27) 
(28) 
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By Theorem [H the term in line (128D is bounded above by 1/ log 2 N, and thus converges 
to 0. For the term in line (f2"7|) . because of the indicator l\v N \<i/2i we use the fact that when 
\ x \ < 1/2, we have \x — log (1 + x)\ < x 2 . Thus this line is bounded above by 



E l/2 



L |yjv|<i/2 



Z(l) 



1/2 N 1-H 



-2,H 



E l/2 



-|VWI<l/2 \ C 2,H 



1/2 N 1-H 



v N \ 



v N 

2 



(29) 
(30) 



The term in line (|29|) converges to by Theorem [3j Finally, by Theorem [T] again, and the 
earlier statement about higher powers of random variables with finite chaos expansions, the 
term in line f)30f) is of order N 2H ~ 2 , and therefore converges to as well. This proves that E2 
converges to 0, finishing the proof of the theorem. 



9.6 Proof of Theorem \7\ 

It is sufficient to prove that 



lim E 

N^oo 



N 



1-H N 



N 1 - 11 } [H N - H)logN 



0. 



We decompose the probability space depending on whether Hn is far or not from its mean. 
For a fixed value e > it is convenient to define the event 



D 



H N > e + 2H - 1 



We have 



E 



N 



E 



+ E 



\.£)C 



N 



N 



N 



l-H N 



1-H 



H 



N 



H) IorN 



N 



N 



l-H 



1-H 



H N -H) logiV 
H n -h) logiV 



+ 



A + B. 



Proof. 



Term A : 

Introduce the notation x = max ( 1 — H, 1 — Hn ) and y = min ( 1 — H, 1 — Hn 



N 



l-H N 



N 



l-H 



^xlogN _ e y\ogN _ e ylogN ( e (x-y) log N _ ^ 



< N y (logN)(x - y)N x ~ y = 2 log NN* 



H -H 



N 



log NN a 



H - H 
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Thus, 



A < E 
= E 



H - H 



N 



log 2 N 



n x-(2-2H) 1dN 2-2H 



H~H N 



' log 2 TV 



Now, choose e € (0, 1 — H). In this case, if lo £ D and x = 1 — H, we get x — (2 — 2i?) = 
—x < — e. On the other hand, for w £ D and x = 1 — Hn we get x — (2 — 2i7) = 
2 — 2iJjy — (2 — 2H) < —e. In conclusion, on D, x — (2 — 2i7) < — e which implies 
immediately 



A < iV~ £ E 



N 



2~2H 



Hn — H 



' log 2 N 



and since the last expectation is bounded 



lim A = 0. 

Term B : 

Now, let uj G D c then # - i^y > 1 - H - e. Since e < 1 - H it implies # > .Hjv- 



Consequently, it is not sufficient to bound 



iV 1 -^ - N 1 - 11 ] above by iV 1 -^. In the 



same fashion we bound 
and | 



H — H 



above by H. Using Holder's inequality with powers \ 



1 D cN 1 ~ Hn 



B < #logiVE 

<H log N[P(D C )] 3/4 ( E 
By Chebyshev's inequality, we have 



1/4 



E 3/4 



p3/4 pc] < 



if -if 



[l-H -£ 



,3/2 



i < cA r- 3 ( 2 - 2H )/ 4 



(31) 



for some constant c depending only if. Dealing with the other term in the upper bound 
for B is a little less obvious. We must return to the definition of if. We have 



Therefore, 



Finally, we get 



1 + V N = N 2 ( H -") = N< H ~") = nK 1 '*) N^- H \ 

E l/4 [^(1-^)41 < N l-H B l/(4) [x + Vn] < 2N l~H 



B < 2Hc(\ogN)N- (1 - H) . 
Finally, B goes to as N — > oo. This finishes the proof of the theorem. 
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